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I The risk minimizing problem 'E[l{{H — min in the Black-Scholes 

OO I framework with correlation is studied. General formulas for the minimal risk function 

and the cost reduction function for the option H depending on multiple underlying 
are derived. The case of a linear and a strictly convex loss function I are examined. 
I Explicit computation for l{x) = x and l{x) = x^ , with p > 1 for digital, quantos, 

. outperformance and spread options are presented. The method is based on the quantile 

• ' hedging approach presented in [1] , [S] and developed for the multidimensional options 
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CN ! 1 Introduction 

O _ 

The paper is devoted to the stochastic control problem arising in the risk analysis of 
financial markets. Let H he a. random variable representing future random payoff which is 
traded on the market. Denote its price determined by the no arbitrage method by p{H). 
If the initial capital x of the writer exceeds p{H) then he is able to hedge H perfectly, 
i.e. he can follow some trading strategy vr such that the wealth process at the final time 
is greater than H, i.e. 

PiX^'"" >H) = 1. 

If X < p{H) then the above equality fails for each vr and as a consequence a shortfall risk 
appears. The aim of the trader is to find a strategy which is optimal in a sense. Let 
I : [0, +oo) — > [0, +cxd) be a loss function which describes the attitude of the trader to 
the hedging losses. The goal is to minimize the shortfall risk defined as 

B[i{{H-x'rn 

This problem was studied with various model settings in many papers. The ones men- 
tioned below do not form a complete list. Existence of the optimal strategy for the case 
when l{x) = x in the context of complete market with the stock prices modeled by the 
diffusion processes was shown in [3] . These results were generalized to incomplete markets 
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in [2] where the existence of solution with the use of dual methods was shown. Existence 
of the optimal strategy in a general semimartingale model was shown in [lOj. More explicit 
results were presented in [S], where the quantile hedging methods, which were introduced 
in [3], enabled to obtain a more precise description of the solution. This paper is con- 
ceptually close to the latter approach and is based on the results obtained in [5j and 
their adaptation for the multidimensional market presented in [1] . We work with the mul- 
tidimensional Black-Scholes model with a correlated Wiener process. A great advantage 
of such a model is its tractability - all parameters can be easily estimated from data, see 
[6], p. 104. It is also complete and one can find explicit formula for the density of the 
martingale measure, see for instance [1]. This enables us to apply the quantile hedging 
methods from [3], [5], [I] which are based on the Neyman- Pearson technique. 

In this paper we study two aspects of the risk minimizing problem for derivatives based 
on multiple underlying. The first is to minimize the risk for a given initial capital x > 0. 
If X > p{H) then the risk can be eliminated by a replicating strategy and thus the minimal 
risk function which will be specified later, equals zero, i.e. ^{{x) = 0. In the opposite 
case ^\{x) is strictly positive and the problem is to find a precise value of ^[{x) and the 
corresponding risk minimizing strategy. The second aspect of the problem is to minimize 
initial costs for the investor who accepts some level of risk v > 0. If u = then the cost 
minimizing function <I>2 equals to the price of H, i.e. $2('^) ~ Pi^) but can be strictly 
smaller if v > 0. The problem is thus to determine ^2('^) ^^"^ ^^"^ '^^^^ minimizing 
strategy. Let us stress the fact that both functions $^,$2 reflect the interplay between 
hedging risk and trading costs and thus serve as important tools for applications. The aim 
of this paper is to present explicit computing methods for the functions ^2 so that to 
be close to the practitioners' needs. 

The idea of the paper is motivated by the formulas presented in Section 6 of [5j which 
concern the call option in a one dimensional Black-Scholes model. The key observation is 
that the solution to the risk minimizing problems can be formulated with the use of two 
real valued deterministic functions. We generalize this concept to the multidimensional 
setting and show that the risk minimizing problems can be solved in the same way as 
well provided regularity of the auxiliary functions. For the linear case we introduce real 
valued functions ^'i,^'2, see formulas (j3.4p . (j3.5p . and show that they are continuous if 
the Wiener process is not degenerate, see Lemma 13.11 and Corollary 13.31 Moreover, for 
many derivatives ^1, ^'2 are strictly monotone, see Lemma [3.H Example (c) and comments 
proceeding formulation of Theorem 13.41 Thus ^'^'^,^'2^ exist and roughly speaking 

^[ = '^lO\if2^7 <l>2 = ^'2 o 

up to the discounting factor and the shift parameter, for precise formulation see Theorem 
13.41 For the general case of a strictly convex loss function / we introduce the functions 
^''i,^'^, see l^m . ([3:23]) . They are regular and Theorem 13.51 , which is actually a refor- 
mulation of Theorem 3.2 in |5J, yields 

^[ = ^[ O {¥2)'^. 

The characterization of the function $2 requires proving of the auxiliary result which is 
formulated as Theorem 13. 6i That result can not be proved with the same method as in 
[5] since the constraints in the associated problem ()3.26p are not linear and this excludes 
possibility of applying the Neyman-Pearson lemma. We present the proof which is based 
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on the Lagrange multipliers. Finally, in Theorem 13.71 we show that 

To summarize, if one is able to find the auxiliary functions ^1, ^2' then he is able 

to find the corresponding functions ^>2- We present explicit computation of the auxil- 
iary functions for several derivatives which are widely traded like digital option, quantos, 
outperformance and spread options. For these examples the auxiliary functions are ex- 
pressed as integrals of the normal densities with appropriate parameters and can effectively 
be applied in practice. 

The paper is organized as follows. In Section [2] we describe the model settings and 
strictly formulate the problem. Section [3] contains the main results which consist of two 
parts concerning a linear and a convex loss function respectively. Section U] is devoted to 
presenting explicit calculations for two dimensional model when l{x) = x and l{x) = ^ 
with p > 1. 



2 Problem formulation 

We work with a multidimensional stock price model with dynamics of d stocks given by a 
standard Black-Scholes model 



dSi = Sliaidt + aidWi), 



1,2,..., d, te[0,T], 



where ai £ R, ai > 0, i = l,2,...,d. Above Wt = {W^ ,Wt^ , ...,Wf),t G [0,T], is a 
sequence of correlated standard Wiener processes. The correlation matrix Q oiW, which 
is assumed to be positive definite, is of the form 



Q 



1 Pl,2 Pl,3 • ■ 
P2,l 1 P2,3 ■ ' 

Pd,l Pd,2 Pd,3 ■ ■ 



Pl,d 
P2,d 



where 



PiJ 



cor 



hj = 1,2, ...,d. 



The process W as above is called a Q- Wiener process. The dynamics the a money market 
account is given by 

dBt = rBtdt, t G [0,T], 

where r stands for a constant interest rate. It is known that such a market is complete 
and that the unique martingale measure P is given by the density 



dP ^ -(Q-M^ 



dP 

with the notation 



Zt := e 



(2.11 



a — rln 



■■=Q-' 



ai —r 

02— r 
o"2 



t€[0,T], 
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for more details see, for instance, [T]. Moreover, 

Wt:=Wt + ^^^t, te[0,T], 
a 

is a Q- Wiener process under P. The dynamics of the prices under the measure P can be 
written as 

dSf = Siirdt + aidWi), i = l,2, d, t G [0,T]. 

The wealth process corresponding to the initial endowment x and the trading strategy vr 
is given by 

1=1 

Each strategy is assumed to be admissible, i.e. X^''^ > for each t G [0, T] almost surely 
and self- financing, i.e. 

d 

dX^'^ = TT^dBt + <dSi , te[0,T]. 

i=l 

A contingent claim is represented by an J^t- measurable random variable H which is 
assumed to be nonnegative, i.e. H >0. As the market is complete, the price of H defined 

by 

p{H) := inf {x : Bvr s.t. > H) = l} 

is given by p{H) = E[e^''^i?], where the expectation is calculated under the measure P. 
The aim of the trader is to minimize the shortfall risk defined by 

where / : [0, +00) — > [0, +00) is a loss function which is assumed to be increasing with 
l{0) = 0. It is clear that if x > p{H) then the risk equals zero for the replicating strategy. 
In the opposite case the risk is strictly positive and the question under consideration is to 
find a strategy such that 

E[/((i/ - ^min. 
We will refer the corresponding function <I>i : [0, +00) — > [0,'E[l{H)]] given by 

<^{{x) := m.mE[l{{H - (2.2) 

as the minimal risk function. The strategy vr such that E[l{{H — Xi^'^)~^)] = ^[{x) will 
be called the risk minimizing strategy for x. If a; > p{H) then <I>'^(a;) = and ^\{x) > 
otherwise. 

We also consider the cost reduction problem. Let v > be a fixed number describing 
the level of shortfall risk accepted by the trader. We are searching for a minimal initial 
cost such that there exists a strategy with the risk not exceeding v, i.e. 

x — >mm; 3 vr s.t. B[l{{H - X^''') + )] < v. 
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The cost reduction function ^2 : [0, +00) — > [0,p{H)] is thus defined by 

^•^(u) := mill {x : 3tt s.t. B[l{{H - < v} . (2.3) 

The strategy vr such that E[Z((X*^^^''''^ — H)^)] < v will be called the cost minimizing 
strategy for v. Notice that ^2(.^) — P{^)- 

In the sequel we examine two cases of the loss functions, i.e. l{x) = x and / a gen- 
eral strictly convex function. The aim is to provide explicit computing methods for the 
functions 

3 Main results 

3.1 Linear loss function 

In this section we examine the case when l[x) = x and denote the corresponding functions 
$2 by $1, $2 respectively. It turns out that the functions ^>i, $2 can be characterized 
in terms of two auxiliary functions 

^i(c):=E(m^J, (3.4) 
^2{c) ■.= BiHUj, (3.5) 

where 

Ac := {Z^^ > c}, c > 0, 



and Zt is given by ()2.ip . 

Let us start with an auxiliary result which establishes regularity properties for the 
functions ^'i, ^'2. 

Lemma 3.1 Let X > 0,Y > be random variables such that EX < +00. Then the 
function g : [0, +00) — )■ [0, +00) given by 

g{c) := E[Xl|y>,|] 

a) is left continuous on (0, +00) with right limits on [0, +cxd), 

b) is right continuous on [0, +00) if the distribution function ofY is continuous, 

c) is strictly decreasing if for any < a < b < +00 holds 

P{X >0,Y e[a,b)) >0. (3.6) 

Proof: The function g is decreasing and thus it has right and left limits. Let us consider 
the auxiliary probability measure P defined by 

'-^ = ^dP 
dP E[X] ' 

which is absolutely continuous wrt. P, i.e. P P. 



a) For any c > we have 

rt ^ ✓ 
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and thus 



I 9{c - ^) - 9{c) 1= E {Xl^^_,^y^^y) = B[X]P (c - 1 < y < c) 0. 
b) For any c > we have 

n{c<y<c + i| = {y = c}, 

and thus 

1)\=e{xi,^^,,^,,,)=b[x]P(c<y<c + 1 ^ 



g{c) - g{c + -) 1= E ( X1^^<^<^+jlj J = E[X]P [c<Y <c+-]^ ¥.[X]P{Y = c) = 0, 
asP <:P and P(y = c) = 0. 

c) Let us notice that p.6p is equivalent to the condition 

3e > s.t. P{X > e, y e [a, 6)) > 0, 
and thus for < a < 6 < +oo we have 

I ^(a) - ai^) I = E (Xl|„<y<b}) = E (Xl|„<y<b}l|j^=o}) + E (Xl|a<y<5}l{x>o}) 
> E (Xl|,<y<b|l|x>,|) > eP(X > e,a < y < 6) > 0. 

□ 

Remark 3.2 Lei ?is notice that the condition ()3.6p implies that Y has a strictly increasing 
distribution function. Indeed, in the opposite case ()3.6p is not satisfied for some < a < 
b < +00. 

Corollary 3.3 // the distribution function of Y is continuous then the function g{c) in 
Lemma \3.1\ is continuous on (0, +oo) and right continuous at 0. 

Examples The condition (|3.6p is satisfied in the following cases. 

a) X > and Y has strictly increasing distribution function. 

b) X, Y are independent, Y has strictly increasing distribution function and P{X > 0) > 
0. 

c) Let {Zi,Z2) be a random vector with nondegenerate normal distribution on a plane. 
Let /, g be functions such that 

/:M2 ^ (0,+oo), 

g : M — > (0, +oo) is strictly monotone. 

Let a, /3, 7, (5 G M be such that the vectors (a, /3), (7, (^) are not parallel, i.e. (a,/3) H 
(7,<5). Then 

X := f{Zi, Z2)l{aZi+i3Z2>k}, y ■= ailZi + 5Z2), 
where k is some constant, satisfy (j3.6p . Indeed, we have 

P{X > 0,y G [a,b)) = p(aZi+l3Z2 > k,g-\a) < 7Z1 + < g-\b) 



for the case when g is strictly increasing. The probabihty above is positive because the 
set 

: {ax + /3y) > k,g~\a) < -fx + 5y < g'\b)^ 

is of positive Lebesgue measure and (^1,^2) has nondegenerate distribution. □ 
Let us notice that due to the fact that Q is nonsingular the random variable 



has a continuous distribution function wrt. P and P. Thus it follows from Corollarv 13.31 
that the functions ^1, ^'2 are continuous. As they are decreasing with images [0,E[ff]], 
[0,e^'^p{H)], the equations 

^i{c)=x, x£[0,E[H]], 

^2{c)=x, xe[0,e'-''piH)], 
have solutions. Moreover, the solution of the first (resp. second) equation is unique if 

P(^H >0,Z^^ £[a,b)^ >0, (3.8) 

resp. 

P(^H >0,Z^^ £[a,b)^ >0. (3.9) 

It follows from Example (c) above that ()3.8p and ()3.9p are satisfied, for instance, when 
d = 2 and 

a) H is a digital option, i.e. H = Kl^gi^yg2^y and (<ti,— (T2) ltQ~^[^~^^] > 

b) H is a quanto domestic option, i.e. H = S^{S^ — K)^ and (cji, 0) H Q^^[^^-p^], 

c) H is a quanto foreign option, i.e. H = (5^ — ■^)^ and (cTi,cr2) Q^^[^^-^^]- 
Below we present the description of the functions <I>i, <I>2- 

Theorem 3.4 a) Let c = c{x) be a solution of the equation 

^f^ic) = e'^^x, X G [0,p{H)). (3.10) 

Then 

['l'i(0)-^'i(c) forx£[0,p{H)), 
\0 forx>p{H). 



^a(x) 



Moreover, the replicating strategy for the payoff Hl^^^^^^ is a risk minimizing strategy 
for X. 

b) Let c = c{v) be a solution of the equation 

^i(c) = ^i(O) - vG[0,E[H]). (3.11) 

Then 

'e-'^^^2(c) forv G[0,B[H]), 



/oru>E[i?]. 



Moreover, the replicating strategy for the payoff HIa^^^^ is a cost minimizing strategy 
for V. 



Proof: For any admissible strategy (x, vr) let us define the success function 

One can check the following identity 

{H - X^''')+ = H - X^'"" AH = H - Hif^^^, 

which implies that 

- = ^H] - nH^xA- (3-12) 

a) In view of ()3.12p the problem ()2.2p of finding $i(x) is equivalent to that of finding the 
strategy vr satisfying 

Y^[HlPx^t,] max. 

If X > p{H) then ipx^-n = 1 for the replicating strategy and <I*i(x) = 0, so consider the case 
< X < p{H). Let us formulate an auxiliary problem of determining ip £ TZ solving 

EilHif]^ max, 

(3.13) 

E[e-''^Hip] < X, 
where 

IZ := {if : < if < 1 and is Tt — measurable}. (3-14) 

It is clear that if (p such that ^[e~'^'^ H0\ = x is a solution of (|3.13p then the replicating 
strategy vf for the payoff H(p is a risk minimizing strategy for x and 

^'i(x) = E[{H - = E[i7] - E[H0\. (3.15) 

Thus now let us focus on determining solution ip of ()3.13p . To this end introduce two 
probability measures Pi , P2 with densities 

dPi H dP2 e-^'^ZxH 



dP E[F]' dP B[e-^^ZTH]' 

Then ()3.13p reads as 

' E-^i [p] — > max, 

(3.16) 

which is a standard problem in the theory of statistical tests. One should try to search 
for the solution in the class of — 1 valued functions of the form Iaci c > 0, where 

UP2 ~ I \dPdP2- I XeIH] ZtH - I I ^ - E[ZtH]I 

For the sake of simplicity we can reparametrize Ac by denoting the constant '^^ above 
just by c. Then Ac is of the form 

Ac :={Zyi >c}. 
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It is known by the Neyman-Pearson lemma that if there exists c = c{x) such that 

E^niAj = P2(^c) = (3.17) 

then the solution of (|3.16p . or equivalently (jS.lSp . is given by = ^A^^^y But let us notice 
that ()3.17p is equivalent to the following 

^2(c) = e'''^x, 

and the existence of the required constant c follows from (|3.1Up . Finally, coming back to 
(|3.15p and using definition of ^'i, we obtain 

$i(x) = B[H] - B[H^] = E[H] - E[H1aJ = ^i(O) - ^i(c). 

b) V > E[if] then the cost minimizing strategy is trivial, i.e. = 0,7r = 0) and thus 
^2{v) = 0. Let us focus on the case when v G [0,E[i/]). In view of (j3.12p the risk 
minimizing strategy is the one which solves the problem 

'B[Hip.,,^]>B[H]-v 

E[e~^'^ Hipx^Tv] — > min. 
We are thus looking for a solution ip £ TZ oi the problem 

'B[Hip] > E[H] -V 

(3.18) 

E[e ^'^Hip] — > min. 

If ()3.18p has a solution satisfying 'E[H(f] = E[ff] — v then the cost minimizing strategy is 
the one which replicates H0 and the cost minimizing function equals 

$2(r) = e-''^E[F(^]. (3.19) 

Let us focus on determining the solution if of (j3.18p . Using notation from the part (a) we 
can reformulate (|3.18p to the form 

(3.20) 



mm . 



It can be shown in the same way as in the proof of Neyman-Pearson lemma that the 
solution should be searched in the — 1 valued functions of the form 1^^; c > 0, where 

IdPi - / XdPdPi- } I ^ ~ cEIZtH]) 



Denoting, for simplicity, the constant - ~ t above by c, we have 

Be = {Zt^ > c}. 
If there exists constant c = c{v) satisfying 



E^MlBj=Pi(i?c) = ^^ (3.21) 
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then (f = is a solution of p.20p or, equivalently, (|3.18p . Let us notice that (|3.2ip can 
be written as 

qi^(c) = ^'i(O) -u 

and existence of the required constant c{v) fohows from ()3.1ip . Coming back to ()3.19p we 
obtain 

□ 

3.2 Convex loss function 

In this section we study the case when / : [0, +00) — > [0, +00) is an increasing, strictly 
convex function such that 1(0) = 0. We assume that / G C^(0, +00) and that /' is strictly 
increasing with ^(0+) = 0, /'(+00) = +00. The inverse of the first derivative will be 
denoted by /, i.e. 

i={ir'- 

The functions $2 can be characterized in terms of the functions 

^[{c):=Bm-^,)H)] (3.22) 

¥^{c):=-E[Hip,]. (3.23) 

where ipc is defined by 

</'c := 1 1 - (^^^^ A 1^ I 1{H>0} , c > 0. (3.24) 

It was shown in [5], Theorem 5.1, that the problem of determining ^\ is equivalent to 
finding the solution ip of the problem 

'E[/((l-v^)i7)] — >min 

(3.25) 

B[e-'^^HyD] < x, 

where 7^ is defined in (fXTil) . Then ¥^{x) = E[/((l - (f)H)] and the risk mmimizmg 
strategy is the one which replicates Hif. Moreover, since the function is continuous 
with the image [0,e^'^p{H)], see the proof of Theorem 5.1 in [5j, it follows that for any 
X S [0, e'^^p{H)] there exists constant c such that ^'2(c) = E[if(/9c] = e'^'^x. Such ipc solves 
the auxiliary problem (j3.25p and thus 

$'i(x) = E[/((l-<^,)i/)], 

and the minimal risk strategy is that replicating the payoff Hipc, see Theorem 3.2 in [5]. 
Thus the results from [5j can be expressed in our notation as follows. 

Theorem 3.5 Let c = c{x) be a solution of the equation 

$^(c) = e'^^x, x€[0,p(H)). 

Then 

$/(^)=/*i(c) forxe[0,p{H)), 
1 for X > p{H). 
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Although Theorem 13.51 is only a reformulation of Theorem 3.2 in [5j, it provides an 
effective method for practical applications if one is able to derive the functions 
concrete derivatives. 

In the sequel we will show that the function <I>2 can be characterized in terms of the 
functions as well. It is easy to show that the cost reduction problem is equivalent 

to that of finding ip £ TZ such that 

Em-ip)H)]<v 

(3.26) 

F,[e~^'^Hip] — > min. 

Let us notice that ()3.26p can not be solved with the same method as ()3.25p . In (j3.25p the 
constraints are linear and thus the solution could be found via Neyman-Pearson approach 
to the variational problem, see the proof of Theorem 5.1 in [5j and p. 210 in [9J. The 
constraints in (|3.26p are no longer linear and the method above fails. Below we present 
the proof based on the Lagrange multipliers. 

Theorem 3.6 The solution of the problem (|3.26p is of the form 
where c is such that E[/((l — <f)H)] = v. 

Proof: First let us notice that if (/? G 7^ is a solution to (|3.26p then necessarily E[/((l — 
(p)H)] = V. Indeed, assume to the contrary that ip is a solution to ()3.26p with E[/((l — 
'p)H)] < V and consider a family of random variables (fa := (f Aa;a £ [0, 1]. Then the 
function a — t- E[/((1 — ipa)H)] is continuously decreasing from E[/(i?)] to 0. Thus there 
exists a G [0, 1] such that E[/((l - (pa)H)] = r. Then ifa < ^ and thus E[Hipa] < E[H(p], 
which is a contradiction. 

Let (f ^ (f he any element of TZ such that E[Z((1 — ip)H)] = v. We need to show that 
E[H(p] < E[Hip]. Let us define ip^ by 

ipe ■■= {l-£)'P + £f, eG[0, 1], 

and the function 

F^{e) := E{HVe) = EiZrHip,). 

We need to show that F^{0) < Fip{l). We will show that has minimum at 0. Let us 
define the auxiliary function 

G^{e) :=Em-^,)H)], 

and notice that due to the convexity of / we have G^{e) < v for each e G [0, 1] 
problem of minimizing F^ on [0, 1] is equivalent to the following 

F^{e) — > min 

< 

e > 0, 
1 - e > 0. 



. Thus the 



(3.27) 
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Both functions F^, G^p are smooth with 

G'^{e)=n'{{^-Ve)H)-{^-v)H], 

and thus the Lagrange function for ()3.27p is of the form 

L{e, Ai, A2, A3) = F^{e) - \i{v - G^{e)) - A2e - A3(l - e). 

As the function F^p is hnear, it attains its minimal value at or 1. We will show that the 
first and the second order differential conditions are satisfied for e = 0. 
The first order conditions are 

Ai, A2, As) = E[Zt(v3 - ^)H] + AiE[r((l - ipe)H) ■ - ip)H] - A2 + A3 = 

(3.28) 

Ai,A2,A3>0, \i{v - G^{e)) = 0, Ase = 0, A3(l-e) = 0. (3.29) 
By the definition of (f we have 

^ = 1- and cZt = /'((I - '^)H) on A 

H 

ip = on A", 

where A := {cZt < l'{H)} and A^ stands for the compliment of A. For e = it follows 
from ()3.29p that A3 = and the equation ()3.28p is of the form 

E[Zt((^ - ^)H1a] + E[Zt{^ - ip)HlA-] + cXiE[Zt{^ - (^)FIa] 

= (1 - cAi)E[Zr((/? - 'f)HlA] + B[Zt(pH1ac] - Xi'E[l'{H)ipHlAc] = A2. (3.30) 
The left side of ()3.30p satisfies the following estimation 

(1 - cAi)E[Zt((/? - (p)H1a] + E[Zt(/?F1ac] - Xi'E[1'{H)lpH1a'^] 

> (1 - cAi)E[Zt(93 - >f)HlA] + B[ZTipHlA-] - XicE[Zt^H1a-] 

> (1 - cAi)E[Zt(93 - <f)HlA + Zt^HIa':]. 

If E[Zt{(P — (f)HlA + Zt^HIa'^] > then we take Ai such that (1 — cAi) > 0, in the 
opposite case, such that (1 — cAi) < 0. In both cases A2 given by ()3.30p is nonnegative. 
The second order condition for e = is 

L'^ie, Ai, A2, A3) = AiE[r'((l - ^)H) • - ipfH^] > 0, 

and thus the solution of (|3.27p is e = 0. 

□ 

The Theorem 13.61 and the definitions of l^^d us to the following result. 

Theorem 3.7 Let c = c{v) be a solution of the equation 

Then 

'e-^^f^(c) forve[0,E[l{H)]), 
forv>E[l{H)]. 
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4 Two dimensional model 



In this section we determine exphcit formulas for the functions "^[^ when l{x) = x and 
l{x) = > 1- In the latter case we use the notation ^'j' = ^5^, ^2 = ^2- examine 
several examples of popular options. 

XP'^ and in view of ()3.24p the following holds 



For l{x) = ^ we have I{x) 



1. 
p 

E 



E 



+ 1e 

P 



H - [cZt] ) 1 



where 



Ac := {cZt < HP-^}, 



(4.31) 
(4.32) 

(4.33) 



and A"^ stands for the compliment of Ac. 

Since our formulas are expressed in terms of integrals of normal densities, at the 
beginning we recall basic properties of the multidimensional normal distribution. They 
can be found in standard textbooks on probability theory or statistics, see for instance [7]. 
A random vector X taking values in has a multidimensional normal distribution if its 
density is of the form 



fx{x) 



(27r)2(detS) 



e 2 



X £ 



(4.34) 



where m S M'^ is a mean of X and S is a symmetric positive definite d x d covariance 
matrix of X. The fact that X has a density (|4.34|) will be denoted by X ~ Nfi(m, S) or 
C{X) = N(i{m, E). If d = 1 then the subscript is omitted and N{m, a) denotes the normal 
distribution with mean m and variance a. If X ~ Nd{m, S) and A is a k x d matrix then. 



AX ~ Nk{Am,AJ:A^ 



in particular if a € M'^ then 



a' X ~ N{a^m,a^ T.a 



(4.35) 



(4.36) 



Let X be a random vector taking values in and fix an integer < k < d. Let us divide 
X into two vectors X^^^ and X^'^^ with lengths k, d — k respectively, i.e. 

X^^^ = {Xi,X2, ...jXk)'^ , X^"^^ = {Xk+i,Xk+2, Xd)'^ . 

Analogously, divide the mean vector m and the covariance matrix S 



m 



so that EX(i) = m(i) , EX(2) = m^^) , CovX^^) = , CovX(^) = S(22) , Cov{X(^) , ) = 
5](i2) ^ Denote by C (X^^) | X^^^ = x^^)) the conditional distribution of X^^) given 

jf(2) ^ ^(2) g ^d-k_ jf 5.(22) -g nonsingular then 











[ m(2) ) ' 


S = 









£ (X« I X(2) = = Ar,(^{i)(a;(2)),s(ii)(^(2))), 



(4.37) 
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where 



= - s(i2)s(22)-^s(2i). (4.38) 

Actually the conditional variance Y,^^^^ (x^^^) does not depend on x^^^ but we keep the nota- 
tion for the sake of consistency. The conditional density will be denoted by fxm\x(^)=x(^^ i^^^^ ) ' 
where x*-^^ G M'^. In particular if {X, Y) is a two dimensional normal vector with parame- 
ters 



then 



where 





( mi \ 




0-11 


0'12 


m = 




S = 








\ m2 J ' 




0'21 


0"22 



C{X\Y = y) = Nimi{y),aiiy)), 



mi{y) :=mi + —{y-m2), o-i{y) := au - — . 

C22 ^22 



(4.39) 



If X is a random vector then its distribution wrt. the measure P will be denoted by C{X) 
and its density by fx- Analogously, fx(i)\x(2)=x('2){x^^^) stands for the conditional density 
with respect to the measure P. 

Below we simplify the multidimensional notation to the case d = 2. The correlation 
matrix is of the form 



Q 



1 P 
p 1 



and thus we have 



1 



-1 p 
P -1 



1 



1 



1 



i-p 
1 



1 



+ 



1 



VO+P v^T+p ' yO^ 

Hence the density of the martingale measure (j2.ip can be written as 



-AiW^-A2W^-BT 



-AiW^-A2W^-BT 



(4.40) 



where 



A, 



A. 



1 



B 



- 1 
1 

1 



ai — r a2 — r 
+ p 

CJi CJ2 

a\ — r 02 — r 



0-2 



1 



+ 



1 



Oil 



0-1 



+ 



1 



1 



a2-r 
02 



+ 



1 



1 



B := B- A 



ai — r 



+ 



A, 



0-1 



0"! 

02 - r 

0-2 



1 



+ 



1 



02 - r 

0-2 
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In the following subsections we will use the universal constants: Ai,A2,B,B defined 
in (|4.4U|) as well as ai,a2,b,ai,a2,b introduced below. 
Fix numbers K > 0, c > 0. One can check the following 

{St >K} = {W^ > ai} = > 5i} , (4.41) 

{St >K} = {W^ > aa} = {w^ > aa} , (4.42) 

{S^ > S^} = {aiW^ - fJaWl >b} = ^aiW^ - aaW^I > ft} , (4.43) 
[Zt^ > c} = {AiW^ + A2W^ > Inc-ST} = ^AiW^ + A2W^ > In c-^t} , (4.44) 



where 



1 ^ ' 1 2^^^ , ._ 1 . I,, 



6 := In (I) + ia2 -a,- \{al - al))T, b := In (|) - i(ai - a?)T. 

In all the formulas appearing in the sequel it is understood that In = — oo and $ stands 
for the distribution function of A^(0, 1). 

4.1 Digital option 

Digital option is a contract with the payoff function of the form 

= K • l|5i>52|, where K > 0. (4.45) 

Let (X, Y), {X, Yjhe random vectors defined by X := aiW^-a2W^, Y := 11^+^2 W|, 

X := piWt — cr2WT, Y := AiWt + ^aW'y- They are normally distributed under P, resp. 

P and their parameters are given by (|4.35|) . 

Linear loss function 

Using ()4.40p and (|4.43p we obtain 

^-1(0) = KE(l|5i >c;2}l^^-i>^p = KP{aiW^ - a2W^ > b, AiW^ + ^2W| > Inc - BT), 
and thus 

/•+00 /•+00 

*i(c)=ic/ / fxx{x,y)dydx. 

Jb Jlnc-BT 



Analogous computation yields 

/■+00 /■+00 

^2{c) = KP{aiW^-a2W^>b,AiW^ + A2W^>lnc-BT) = K / fx,Y{x,y)dydx. 

Jb Jlnc-BT 
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Power loss function 

In view of and ([OCT]) we have 



A, := {cZt < H^-^} = {cZt < K^-^^^^i.^^w^^.y} = {a^W^ - a^W^ > b, cZt < K^-'} 
= {aiW^ - daWl > b, AiW^ + AiW^ > In i —— J - BT} (4.46) 



{aiW^ - £72 W| > b, AiW^ + A2W^ > In 



KP- 



BT}, 



and thus 



(4.47) 



r,lAj-C—^B[Z^-'lAj. 



In view of (j4.46p and (j4.47p we have 

RP f f KP^^\ \ 1 p ^- 

$P(c) = —P / - (y2W^ > b, AiW^ + AiW^ < In ( — ^ ) - BTj + -c~E[Zj^-' IaJ 



P Jb J-oo 



p /■+00 /.+00 

fx.Y{x,y)dydx + -c~ 



p(y+BT) 

e fx,Y[x,y)dydx, 



and 



K 



oo /" + 00 



fxY(^,y)dydx - Cf-i 



+ 00 /■+00 



y+BT _ 



4.2 Quantos 

4.2.1 Quanto domestic 

The contingent claim is of the form 

H = SUSt - K)+, K >0. 

Linear loss function 

Using (|4.4ip we obtain 



(4.48) 



^'i(c) = B[SUs't - i^)+l|^^i>,}] = eIsUs't - ^)l{z->c} \S't>k] P(4 > K) 



E 



SUSt - K)l{A^W^+A,W}>lnc-BT} I > ^l] ^(W^ > ^l) 



+00 



E 



0.1 



J ai 



+00 



fwi{x)dx 



^) /inc-sT-Ai. ^^'^ fw}\w^=xiy)dyfw^{x)dx, 
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and 



E 



S2g{r-|^i)T+a2VK|(5lg(r-l<72)r+aiH7^ - K)l 



+00 



E 



ai 



J a\ J J- 



Power loss function 

The set ()i33]) is of the form 

Ac := {cZt < HP-^} = < 



ce 



-AiW^-A2W}-BT\ p-1 



< (5-]. - KY 



cp- 



02 



For simphcity we assume that + cr2 > 0. In the opposite case one has to modify the 
form of the set Ac and thus also the integration hmits in the formulas below. We obtain 

Ac = {W^ > w{W^),W^ > ai} = jw^l > w{W^),W^ > di} , 



where 



■w{x) :- 



■w{x) :- 



cP- 



-^x + In 
1 



1 

CP-l 



In view of this above, (j4.3ip . (|4.32|) and using conditional densities we obtain 



V 

+00 f+00 



00 /"+00 



ai 



+ 



ai J w{x) 
p _ BTp 



00 r+oo 



T2X _ 



Wj,\W^=x 



{y)fw:^{x)dydx 



M/P(c) = 52e('-5-i)^ / 

J a 



a\ Jw(x) 
+00 f+OO 



1 BT 

cp-^e 



ai J w{x) 
+00 /"+00 



ai J w(x) 



Kf fw^\w^=x{y)fw^{x)dydx 
e vP-i p-^yj f^2\w^^x{y)fw^{x)dydx, 
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4.2.2 Quanto foreign 

The payoff is of the form 

K 



H = [S\ - — ] , K>0. 



q2 



Linear loss function 

First let us notice that 



- ^ > o| = {aiW^ + a2W^ >d} = {cxiVf^ + (^2^1 > d} , (4.49) 



where 



d ■■= In - («i + «2 - ^(^7? + al)^ T, d := In - (^2r - ^(^? + ^l^j T. 

(4.50) 

We have 



^i(c) = E 



E 



9^ ^Vi 



- ^) 1 , lne-BT-A,wl , I ^1 + (72^1 > d] P{aiW^ + CJ2W^ > d) . 

bj,J {wj,> J-} J 



Denoting Z := aiW^+a2W^ and taking into account conditional distribution £(Ty^, \ Z) 
we obtain 



-OO /'+00 /' + 00 

*i(c) ^J^ J ^ /„^_,,_,^J gcle("^-^-^)^+-^- - i^5o^e(-"^+^-^^)^-'^^^)/(^.,^.)l^^,(x,y)dydx/z(z)dz. 



>l2 

Using the same argument under the measure P with Z := aiW^ + a2W^ yields 

*2(c) = e[ - ^) l^^^^,,,^_sT^^^^^^ I aiW?^ + aaW^I > d\p{aiW^ + a2W} > d) 

I \ brp J {VVj,> } J 

+ 00 /' + 00 /' + 00 



In c-BT-A-i 



Power loss function 

Using ()4.49p one can check the following 

A.:.{.Z..((s|--)7-.S|-|>0} 



+ -£72)1^1 > v{aiW^ + a2W^),aiW^ + a2W^ > d] (4.51) 
Ip — 1 p — 1 ) 

{^^W^ + {^^-a2)W^ > v{aiW^ + a2Wh,aiWl + a2W^ > d\, (4.52) 
Ln— l^p— 1' J 
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where d, d are given by (j4.50p and 



, f 5l52e(-l+-2-i(^?+<xi))T+x_^■ 
^' (X) = m < - 



CP-IDqC^ ^22 p-lJ 



vix) = m < = 

CP-ISqC^ 2 2 p-l-l 

To calculate ^'^^,^'2 use conditional distributions C{X \ Y), C{X \ Y), where X := 
^Wt+{^-<'2)WI Y := aiW^ + a2Wl Y := -.1^^^^ + (-^ - ^2) y - + 

(72 W^. Denote by fci, /c2, ^3, ^4 constants satisfying = /ciX + /c2^5 ~ + k^Y, 
= kiX + k2Y, = hX + k^Y. Then we have 



1 /-oo /.ti(i/) / 



^lg(Q2-5Cr|)T+<Ti(fc3X + fe4t;) 



.fx\Y=yix)fY{y)dxdy 



1 P pBT 



00 /'+00 



/'+00 p-\-oo I TC \ 

*2(c) - / / S'ie('-5-?)7'+-i(fei-+fe2y) ^ ^ — U.y(x)fy{y)dxdy 



1 BT 
Qp-1 ^ p-1 



00 /'+00 



4.3 Outperformance option 

The problem is studied for 

H = (max{5^, S^} - K)^ , K > 0. 

Linear loss function 

By (031]), (I02]l and (lOall we get 



+ E 
= E 
+ E 



(S^ - K)l^^-i^^^ \S^>K,S^> > K, > S^) 

{S^ - K)l^^-i^^y I > ai,aiW^ - a2W^ > b P{W^ > ai,aiW^ - a2W^ > b) 
{Si - I W| > a2,aiWl - a2W^ < b P{W^ > a2,(7iWl - a2W^ < b) 



+00 f+CO 



+ 



0.1 

/•+00 fb 
1 02 



(Soe^"! a'^JT+aix -K)l.j^ a^^^^^^_^^^f^i^^^y^i_^^y^2{x,z)dzdx 

(72 ~ 



2la2-\ul)T+cj2y 



K)l 



{Ai 



-A2y>lnc-BT}fw},'^iW^-a2W^iy^z)dzdy, 
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+ E \{S^ - K)l^z-^^c] \S'^>K,S^< S^] P{Sl >K,S^< S^) 

= E {S^- K)1^2-^^c} I > ai,aiW^ - a2W^ > b P{W^ > di,aiW^ - a2W^ > b) 
+ E {S^- K)1^2-^^c} I > a2,c7iW^ - a2W^ < b P{W^ > 03, aiW^ - a2W^ < b) 



+00 /"+00 



K)l 



+ 



ai 
+00 

0.2 



{Aix+A2^^^^>lnc-BT}-l'w^,aiW^-a2Wr 



. (x, z)dzdx 



"2 



'?:Ar-Ul)T+a2y 



K)l 



fw^,a,w^-a2W^(y^^)dzdy. 



-00 ' "^1 

Power loss function 

Taking into account ()4.4ip . ()4.42p . ()4.43p we can write 
= {c^T < (^T V - 5^ V 5| - > 0} 

= {cZt < (^T - Kf-\S^ > K,S^> S^} U {cZt < {S^ - Kf-\S^ > K,S^< S^}. 
We consider the case when Ai > 0, > 0: 

Ac = {W} > - (AiW^ + BT + ln (-{S^ - Kf-^)) ) > ai,aiW^ - a2W^ > b} 



U{W^ > -(^A2W^ + BT + ln (^-{S^ - K)P-^)jj > a2,aiW^-a2W^ < b} 
> vi{W^), W^>ai,W^< 

U <! > V2iW^),W} > 02, < 







02 


^} 


a2Wl - 






-} 






0-2 




(72VF|- 


b\ 





(4.53) 



where 



V2{x) = -i- (^A2X + BT + ln (i(52g(a2-iai)T+a2X _ ^ 



vi{x) 

V2{x) 



(^A^x + BT + In (i(5ie(^-^'^?)^+'^i^ - i^f)) , 
(^A2X + BT + In (i(52e('^-^'^')^+'^2^ - i^f . 



Using the representation (|4.53p and accepting the convention that the integral over the 
empty set is zero, we obtain 
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+00 /•vi{x)A 



= / / 

J ai J—oo 



+ 



+ 



<Tix—b 
CT2 

ai Jvi{x) 



Kffw^\w^=xi.y)dyfw^{x)dx 



+00 /•V2{x)A 



a2 •! —00 



(723; — b 



KTfw^\w}=xiy)dyfw}ix)dx 



, CT'jX — b 

+ -cp-^e e p-i p-^yfw^\w^=Ay)dyfw^ix)dx, 

P J ai Jvi{x') 



CT] x — b 
f' + OO p— 

(-) = / / (^o^e(-.-?)^+-^ - K)f^i^^,^^{y)dyf^,{x)dx 

J ai J vi(x) ^ ^ ^ 



(71 x — b 



cp-^e v-^-' j e v-^"" p-^yf^i^^,^^{y)dyf^,{x)dx 



CT2 



(T^yX — b 

r+oo f— 

+ / / (5o^e(-.^^)^+-^ - K)f^,^^i^^{y)dyf^Ax)dx 



(j-^^x — b 

-1 / / /-> T) — 1 n — 1 ^ 



2)T+a2X 



'5,2 Jv2{x) 

4.4 Spread option 

The payoff is of the form 

H = {S\. - Si - K)'^ , K>0. 

One can check the following 

{S^ >Sl + K} = {W^ > d{W^)} = {W^ > d{W^)}, 

where 

d{y) := — In , d{y) := — In ' 

^1 \ S'Ae("i"2'^i)^ / CTi 



Linear loss function 

We have 



r+oo /■+(X) 

-00 Jd{y) 



+00 



/-too p-\-oo . . 
-00 Jd(v) 



and 



/+00 
E[(5^ - 5| - I W^l = y]f^Jy)dy 
-00 



-OO Z' + OC 



-00 Jd{y) 
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Power loss function 

We have 

:= {cZt < (S't -S^- K)p-'^,S^ - - K > 0} 

= {w:^ e Aiw})} = {w:^ e AiW^)}, (4.54) 
where 

A{y) := {x : c^e'^^'^y-^^ < 3^6^''^-^^'')^+'''^ - 52g(a2-|<xi)T+<x2j/ -K,x> d{y)}, 
A{y) := {x : c^e-^^'^^'^^ < 5(5e('-5-?)^+'^i^ - 5^e("-5-i)r+-2S/ - i^, x > d{y)}. 
Let us notice that the set A'^Ci {H > 0} is of the form 

A^^r\{H >0} = {h4 e 13{W^)}, (4.55) 

where 

B(y) := {x : ci^e'F^'^"^^"^^ > ^o^e^^i-^'^i^^+'^i^ - 5^e("2-|<x2)T+<x2y _k,x> d{y)}. 
Taking into account ()4.54p and ()4.55p we obtain 

" J — oo JBiv) 

^''~^^) fw^\w}=yi^)dxfw}{y)dy, 



1 ^ _PST r+f^ /■ / pAi^ pA2. 

H — cp-ie p-i / / (e 

^3 i-oo iyl(y) ^ 



^lic) = j^^^^ (5oie(-l'^?)^+-- - Sle^^-l"l)^^^-y - K)f^,^^^^^{x)dxf^Ay)dy 



1 BT 



+ cf-ie p-i / /_ le f-i"' ''-'^)fwl\wl=J^)dxfw2{y)dy. 

-oo ^ 



' Tl "T" 
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